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Abstract: We summarise a Dyson-Schwinger-equation-based calculation of an extensive range of 
light- and heavy-meson observables, characterised by heavy-meson leptonic decays, semileptonic 
heavy-to-heavy and heavy-to-light transitions - B — > D* , D, p, tt; D — > K* , K, tt, radiative and strong 
decays - B* s ^ — ► -B( s )7; -Df% — * D^ s ^, D?t, and the rare B — > K*^ flavour-changing neutral-current 
process. In the calculation the heavy-quark mass functions are approximated by constants, interpreted 
as their constituent-mass: M c = 1.32 GeV and Mb = 4.65 GcV. 
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1. Introduction 

The Dyson-Schwinger equations (DSEs) [1] pro- 
vide a nonperturbative, Poincare-covariant, field 
theoretical approach to the calculation of hadronic 
matrix elements, and they have been widely ap- 
plied to the phenomena of continuum strong QCD 
[2] . Many applications have focused on nonhadro- 
nic electroweak interactions because the electro- 
weak probe is well understood and the interac- 
tions therefore explore the structure of the hadro- 
nic target. These are just the phenomena of in- 
terest to this community. 

In Refs. [3-6] both light- and heavy-mesons 
are represented as bound states of a dressed-quark 
and -antiquark, with the quarks' dressing des- 
cribed by a DSE: the QCD gap equation. The 
general form of the solution of this equation is 

S f (p) = Z f (p 2 )/[i 7 .p + M f (p 2 )], (1.1) 

/ = u, d, s, c, b is the flavour label, and extensive 
studies have revealed [5,6] that while the mass 
function of a light-quark: M Ui d, s , is a rapidly- 
varying function, that of an heavy-quark is al- 
most momentum-independent: see Fig. |l|. 
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Figure 1: Quark mass function obtained in Ref. [7] 
as a solution of the QCD gap equation with current- 
quark masses, fixed at a renormalisation point £ = 



19 GeV: ml 



3.7 MeV, mi 



82MeV, m% = 



0.58 GeV and m£ = 3.8 GeV. The indicated solutions 
of M 2 (p 2 ) = p 2 define the Euclidean constituent- 
quark mass, Mf, which takes the values: Mf — 
0.56 GeV, Mf = 0.70 GeV, Mf = 1.3 GeV, Mf = 
4.6 GeV. (Figure adapted from Ref. [6].) 



The behaviour of M Cj b (and also that of Z c ^, 
which is not illustrated here) suggests that the 
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heavy-quark propagator can be well-approxima- 
ted by: 

S Q {p) = l/[i 1 -p + M Q ], Q = c,b, (1.2) 

where Mq is a constituent-heavy-quark mass pa- 
rameter. A good description of observable phe- 
nomena requires Mq » Mq, see below. 

The dressed-quark propagator can also be 
written in the form 

s f (p) = -i7-p4ip 2 ) + 4(p 2 )- (i-3) 

and the contrasting, significant momentum-de- 
pendence of the light-quark propagator is effica- 
ciously represented via the algebraic parametri- 
sation introduced in Ref. [8] 

a f s (x) = 2m f T{2{x + rh))) (1.4) 
+ F{b lX )F(b 3 x) (b{ + b f 2 F{exj) , 

^—[l-^filmf))], (1.5) 

f = u,s (isospin symmetry is assumed), T(y) = 
(1 — e~ y )/y, x — p 2 /X 2 ; fhf — m//A; and 

4(x) = \a f s (p 2 ), 4(x) = A 2 4(/), (1.6) 

with A = 0.566 GeV a mass scale. This alge- 
braic form combines the effects of confinement 
and DCSB with free-particle behaviour at large 
spacelike p 2 : <Jv{p 2 ) ~ 1/p 2 and as{p 2 ) ~ m/p 2 . 
(With S(p) an entire function a sufficient con- 
dition for confinement is satisfied; i.e., the ab- 
sence of a Lehmann representation for coloured 
Schwinger functions [2]. e — 10~ 4 is introduced 
in Eq. ( |l.4| ) only to decouple the large- and in- 
termediate-p 2 domains: it is not a fitting param- 
eter.) 

Bound states in quantum field theory are de- 
scribed by Bcthe-Salpeter amplitudes whose mo- 
mentum-dependence implements the necessary re- 
striction of the relative momentum of the con- 
stituents in a bound state. This means that a fur- 



ther, commonly-used approximation to Eq. (1.2) 



S Q (k + P) = ±- 1 



Z7 • v 



2 k-v-E 



O 



H 



Mo 



(1.7) 

where P M =: m H v^ =: (Mq + E H )v^, m H is 
the hadron's mass and k is the momentum of 



the lighter constituent, can only be reliable if 
both the momentum-space width of the Bcthc- 
Salpeter amplitude, lob, and the binding energy, 
Eh, are significantly less than Mq. 

2. Heavy-Quark Limit 

The DSE framework reproduces all the acknowl- 
edged consequences of heavy-quark symmetry. In 
addition, one obtains explicit expressions for the 
physical observables in the heavy quark limit; 
e.g., the leptonic decay constants is given by 



fp = fv 



du 



-E H ) 



V 





(2.1) 

duip 2 H (z) \a f s (z) + y/uo-{,{z)} , 



x ip H (z) <j o- J s (z) + -\/uo- J v {z) 



N_c_ 

4tt 2 



where z = u — 2Eh\/u, f labels the meson's 
lighter quark and ipn(z) is the scalar function 
characterising the dominant Dirac-covariant in 
the heavy-meson's Bethe-Salpeter amplitude; e.g, 
the 75-term for the pseudoscalar and 7 M -term for 
the vector meson. 

As another example, the semileptonic heavy- 
to-heavy pseudoscalar transition form factors (Pi 
— > Pitv) acquire a particularly simple form in the 
heavy-quark symmetry limit: 



f±(t) 



mp 2 ± JTOpj 



£/0), 



(2.2) 



£/0) = K f 



2^/mp 2 m Pl 

1 oc 

N f T f 

' ' du(p 2 H (zw) (2.3) 



4tt 2 J W 


a f s (z w ) 



£4m 



with W 



1 + 2r(l - t)(w - 1), z w 



2E Hy /u/W and 



m. 



m 



w = 



P-2 



2mp 1 mp 2 



= -v Pl - Vp 2 



(2.4) 



The canonical normalisation of the Bethe-Salpeter 
amplitude automatically ensures that 



1) 



1 



(2.5) 
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and it follows [3] from Eq. (2.3) that 



dw 



1 

> - 

- 3 



(2.6) 



Similar analysis for the heavy-to-heavy tran- 
sitions with vector mesons in the final state and 
for heavy-to-light transitions yields relations be- 
tween the form factors that coincide with those 
observed in Ref. [9]; i.e., in the heavy-quark limit 
these form factors too are expressible solely in 
terms of 

3. Results 

This phenomcnological application of DSE meth- 
ods to the calculation of heavy-quark observables 
is founded on a large body of work that has fo- 
cused on light-meson physics. References [2,10] 
provide an overview. Herein we summarise re- 
sults for an extensive but not exhaustive body of 
observables: heavy-meson leptonic decays, semilep- 
tonic heavy-to-heavy and heavy-to-light transi- 
tions - B — * D* , D, p, it; D — > K*, K, tt, radia- 
tive and strong decays - B* s ^ — > £?( s )7; D* s -, — ► 
Dr s yy, Dw, and the rare B — > K*"f flavour-chan- 
ging neutral-current process. 

In Ref. [6] an algebraic characterisation of 
the dressed-quark propagators and bound state 
Bethe-Salpeter amplitudes employing ten param- 
eters, plus the four quark masses, was used in a 
X 2 -fit to iV bs = 42 heavy- and light-meson ob- 
servables. That yielded: x 2 /d.o.f = 1.75 and 
X 2 /N bs — 1.17, and the quality of the fit is il- 
lustrated in Tables [l] and ||. Using an approxi- 
mating algebraic representation of the functions 
involved materially simplifies the calculation of 
observables, bypassing the repeated solving of 
nonlinear, coupled integral equations. It is a use- 
ful but not necessary artefice, as Refs. [7,11-13] 
make plain. 

The fitting yielded dressed-quark-propagator 
parameter values 



m f 



b{ b{ 



0.00948 2.94 0.733 
0.210 3.18 0.858 



(3.1) 



with 6q- s = 0.131, 0.105, by = 0.185, which 
were not varied i.e., retain the values fixed in 



previous studies of light-meson observables [14]. 
The dimensionless u, s current-quark masses in 
Eq. ( |3.l|) correspond to m u — 5.4 MeV and m s = 
119 MeV, and these algebraic propagators yield 
Euclidean constituent-light-quark masses: M,f = 
0.36 GeV and Mf = 0.49 GeV. 

The fitted constituent-heavy-quark mass pa- 
rameters are 

M c = 1.32 GeV and M b = 4.65 GeV , (3.2) 

consistent with the estimates reported in Ref. [15] 
and hence the heavy-meson binding energy is 
large: 



E D := rn D - M c = 0.67 GeV , 
E B := m B - M b = 0.70 GeV . 



(3.3) 



These values yield E D /M C = 0.51 and E B /M b = 
0.15, and provide an indication that while an 



heavy-quark expansion, Eq. (1.7), will be accu- 
rate for the 6-quark it will provide a poor approx- 
imation for the c-quark. The constituent-heavy- 
quark-masses in Eq. ( |3.2[ ), obtained in a Poincare 
covariant approach [6], are, respectively, ~ 25% 
and ~ 10% smaller than the values used in non- 
relativistic models. 

Reference [6] represented the dominant scalar 
function in the light-pseudoscalar-meson Bethe- 
Salpeter amplitude as 

£ P (k 2 ) = -j-B P (k 2 ), P = tt 7 K; (3.4) 
Jp 



constructed via Bp := B u \ b u^ b p using Eq. (]L 

with, e.g., = f n /y/2. This Ansatz follows 
from the constraints imposed by the axial- vector 
Ward-Takahashi identity and the fit yielded 



bl = 0.204 , 



b% = 0.319. 



(3.5) 



The exploration of light-vector-meson prop- 
erties is less extensive than that of light-pseudo- 
scalar-mesons, and this is not peculiar to DSE 
analyses. Hence, following, e.g., Refs. [16], the 
algebraic characterisation of Ref. [6] used 



^ 2 ) = i/(i + fc 4 A4 ); 



(3.6) 



i.e., a one-parameter form to describe the domi- 
nant scalar function in the vector Bethe-Salpeter 
amplitude. This is merely a simple, efficacious 
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u 

fx 

(uu) 1 / 3 

f P 
r 

1 pTTTT 

Id 
fs 



Obs. 

0.131 
0.160 
0.241 
0.245 
0.216 
0.151 

0.200 ± 0.030 
0.170 ± 0.035 



Calc. 
0.146 
0.178 
0.220 
0.255 
0.163 
0.118 
0.213 
0.182 



m K 
(ss 



,1/3 
1/3 



fx* 

r 



K'(Ktt) 
fD. 

gBK'fMb 



Obs. 

0.138 
0.496 
0.227 
0.287 
0.244 
0.051 

0.251 ± 0.030 
2.03 ± 0.62 



Calc. 
0.130 
0.449 
0.199 
0.296 
0.253 
0.052 
0.234 
2.86 



Table 1: The 16 dimension-GeV quantities used in x 2 -fitting the model parameters. The values in the "Obs." 
column are taken from Refs. [7,15,17]. (Table adapted from Ref. [6].) 





Obs. 


Calc. 




Obs. 




Calc. 


/f^(0) 


0.73 


0.58 




0.44 ± 


0.004 


0.44 


Fir (3.3 GcV 2 ) 


0.097 ± 0.019 


0.077 


B(B -> D*) 


0.0453 


± 0.0032 


0.052 


P 2 


1.53 ± 0.36 


1.84 




1.25 ± 


0.26 


0.94 




0.86 ± 0.03 


0.84 


AB^D* 


0.19 ± 


0.031 


0.24 


£(1-2) 


0.75 ± 0.05 


0.72 


B(B -> tt) 


(1.8 ± 


0.6) x 10-4 


2.2 


£(1.3) 


0.66 ± 0.06 


0.63 


fB^TT 

■>+ (14.9 GeV 2 ) 


0.82 ± 


0.17 


0.82 


£(1-4) 


0.59 ± 0.07 


0.56 


J-B^TT 

■>+ (17.9 GeV 2 ) 


1.19 ± 


0.28 


1.00 


£(1.5) 


0.53 ± 0.08 


0.50 


fB — >7T 

•'+ (20.9 GeV 2 ) 


1.89 ± 


0.53 


1.28 


B(B -> D) 


0.020 ± 0.007 


0.013 


B(B -> p) 


(2.5 ± 


0.9) xlo -4 


4.8 


B(D -> K*) 


0.047 ± 0.004 


0.049 


jf^(0) 


0.73 




0.61 


V(0) D ^ K ' 

Mo) 


1.89 ± 0.25 


1.74 




0.73 




0.67 


r T 


1.23 ± 0.13 


1.17 


QB'Bit 


23.0 ± 


5.0 


23.2 


A 2 (0) D ^ K 
Ax(0) 


0.73 ± 0.15 


0.87 


9D*Dir 


10.0 ± 


1.3 


11.0 



Table 2: The 26 dimensionless quantities used in x 2 -htting the model parameters. The values in the "Obs." 
column are taken from Refs. [15,18-22]. The light-meson electromagnetic form factors are calculated in impulse 
approximation [8,14,23] and £(w) is obtained from f+^ D (t) via Eq. ([2.2;). (Table adapted from Ref. [6].) 



Ansatz, which might now be improved by build- 
ing upon Refs. [11,12]. No reliable Bethe-Salpeter 
equation studies exist for heavy-mesons (an effi- 
cacious Bethe-Salpeter kernel is still being sought) 
and therefore the following single-parameter form 
was employed to represent the amplitudes: 



<p H (k 2 ) = exp(-fc 2 /w|f) • 



The x 2- fit produced: 



(3.7) 



(3.8) 



The ordering in magnitude is qualitatively un- 
derstandable: the heavier the meson the smaller 
the spacetime volume occupied. In addition, the 









p 


0.515 


D 


1.81 


K* 


0.817 


B 


1.81 



result: lud = lob, which means that the Comp- 
ton wavelength of the c-quark is greater than the 
length-scale characterising the bound state's ex- 



tent, emphasises that Eq. (1.7) must provide a 
poor approximation for the c-quark. 

With the model's parameters fixed, it is pos- 
sible to calculate a wide range of other light- 
and heavy-meson observables. Some of the re- 
sults are summarised in Tables while Figs. || 
and H depict the calculated ^-dependence of the 
semileptonic transition form factors that are the 
hadronic manifestation of the b — > c, b — > u, 
c — > s and c — > d transitions. The form factors 
can be approximated by the monopole 



h{t) = h(0)/[l-t/h 1 ], 



(3.9) 
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Figure 2: Left panel: calculated semileptonic B —* D and B — > D* form factors. Right panel: the 
semileptonic B — * n and B —> p form factors with, for comparison, data from a lattice simulation [21] and a 
vector dominance, monopole model: /+ Mr (i) = 0.46/(1 — t/m%»), tub* = 5.325 GeV, the light, short-dashed 
line. Monopole fits to the model's results are given in Eqs. (3.9) and (3.10). (Figure adapted from Ref. [6].) 
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Figure 3: Calculated semileptonic D — > K and D — > if* (left panel), D —> n and D — > p (right panel) 
form factors. Monopole fits to the calculated results are given in Eqs. (3.9) and (3.10). (Figure adapted from 
Ref. [6].) 



with h(0) given in Tables [| and ||, and hi, in 
GeV 2 , listed in Eq. ( ^To| ). 







B - 




(4.63) 2 


S - 




(5.58) 2 


L> - 


-> if, if* 


(2.31) 2 


L> - 


-> 7T,p 


(2.25) 2 



/if 1 



hX 



(5.73) 2 (4.64)2 (4.61) 2 
-(21.5) 2 (6.94)2 (7 06)2 
(6.70)2 (3 09)2 (2 78)2 
(7.06)2 (3 17)2 (2 80)2 
(3.10) 

With these calculated results it is possible to 
check the fidelity of heavy-quark symmetry lim- 
its. The universal function characterising semilep- 
tonic transitions in the heavy-quark symmetry 
limit, £(u>), can be obtained most reliably from 



B — > D, D* transitions, if it can be obtained at 
all. Using Eq. (|J) to extract it from ff^ D {t) 
one obtains 



(3.11) 



£'+(!) = 1.08, 



which is a measurable deviation from Eq. ( |2.5| ). 
The calculated form of £*+ (w) (0) is depicted 
in Fig. |] and compared with two experimental 
fits [20]: 



1-p 2 (w-1), 
2 - 0.91 ±0.15 ±0.16, 



P 



i{ w ) = 



exp 



(l-2p 2 )' 



(3.12) 



(3.13) 
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accessible heavy- meson observables [6] . That ex- 
tension explored the fidelity of a simple approx- 



Figure 4: Calculated form of compared with 
experimental analyses. Model: solid line. Exper- 
iment: short-dashed line, linear fit from Ref. [20], 
Eq. ( |3. 1 3| ); long-dashed line, nonlinear fit from 
Ref. [20], Eq. ( |3.14| ). The two light, dotted lines are 
this nonlinear fit evaluated with the extreme values 
of p 2 : upper line, p 2 = 1.17 and lower line, p 2 = 1.89; 
data points, Ref. [25]. (Figure adapted from Ref. [6].) 



p 2 = 1.53 ±0.36 ±0.14. 

The calculated result is well approximated by 
1 



£ /+ (w) 



imation, Eq. (L2), to the dressed-heavy-quark 
propagator, and yields a unified and uniformly 
accurate description of an extensive body of light- 
and heavy-meson observables. Algebraic analysis 
proves [6] that in the heavy-quark limit pseu- 
doscalar meson masses grow linearly with the 
mass of their heaviest constituent; i.e., mp oc 
riiQ, while the numerical results indicate that 
corrections to the heavy-quark symmetry limit 
of < 30% are encountered in b — *■ c transitions 
and that these corrections can be as large as a 
factor of two in c — > d transitions. 

The calculation of the semileptonic transi- 
tion form factors for B- and D-mesons on their 
entire kinematic domain and with the light-quark 
sector well constrained is potentially useful in the 
experimental extraction of the CKM matrix el- 
ements V c b, V u b- That is also true of the calcu- 
lation of the leptonic decay constants; e.g., ac- 
curate knowledge of Jb can assist in the deter- 
mination of Vtd- The calculations show that the 
leptonic decay constants for D /-mesons do not lie 
on the heavy-quark l/-\/mQ-trajectory, and pro- 



l + pj + (w-1) 



Pf + — 1-98 . (3.14) vide an estimate of the total width of the 



The semileptonic B — > D* transition can 
also be used to extract £(w). That yields [6] 



C 1 (1) = 0.987 , C 2 (1) = 1-03 , C (1) = 1.30 , 

(3.15) 



and an w-dependence well-described by Eq. ( 3.14 ) 
but with 

p% = 1.79 , p\ 2 = 1.99 , p 2 v = 2.02 . (3.16) 

These and other results in Ref. [6] furnish a mea- 
sure of the degree to which heavy-quark sym- 
metry is respected in b — > c processes. Com- 
bining them it is clear that even in this case, 
which is the nearest contemporary realisation of 
the heavy-quark symmetry limit, corrections of 
< 30% must be expected. 

4. Epilogue 

Herein we have summarised a direct extension 
of DSE-based phenomenology to experimentally 



and _D*°-mesons, for which currently there are 
only experimental upper-bounds. 

The model we have summarised employs sim- 
ple parametrisations for the dressed-quark prop- 
agators and meson Bethe-Salpeter amplitudes. 
Its efficacy supports the interpretation that heavy- 
and light-mesons are both simply finite-size bound 
states of dressed-quarks and -antiquarks; i.e., they 
are not qualitatively different. Furthermore, this 
efficacy demonstrates that the qualitative elements 
of the Poincare-covariant DSE- framework are rich 
enough to account for the gamut of strong inter- 
action phenomena; i.e. , as observed too in Ref. [26] , 
no essential element is inherently lacking. 

Naturally the model can be improved; e.g., 
via a wider study of light-vector-meson observ- 
ables, so as to more tightly constrain their prop- 
erties, perhaps using direct Bethe-Salpeter equa- 
tion studies like Refs. [11,12] as a foundation 
for improved models of the vector meson Bethe- 
Salpeter amplitudes. Along this path, a more sig- 
nificant extension is the development of a Ward- 
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Obs. 


Calc. 




Obs. 


Calc. 




0.472 ± 0.038 


0.46 


- f 2 r 2 
J is. i\ 


(0.19 ± 0.05) 2 


(0.10) 2 


9 p7T7T 


6.05 ± 0.02 


5.27 


T D , (MeV) 


< 2.1 


0.020 




6.41 ± 0.06 


5.96 


Tn,+ (keV) 


< 131 


37.9 


Qp 


5.03 ± 0.012 


5.27 


T D , D 7 (MeV) 


< 1.9 


0.001 


f D - (GeV) 




0.290 


r B .+ B +^ (keV) 




0.030 


f D » (GeV) 




0.298 


r B »o B o^, (keV) 




0.015 


fB s (GeV) 


0.195 ± 0.035 


0.194 


T B , Bsl (keV) 




0.011 


/b. (GeV) 




0.200 


B(£>*+^ L>+7T°) 


0.306 ± 0.025 


0.316 


f B . (GeV) 




0.209 


B(L>*+^ D°tt+) 


0.683 ± 0.014 


0.683 


W/o 


1.10 ± 0.06 


1.10 


B(L>*+^ L»+ 7 ) 


U.U11 -0.007 


0.001 


/s s / /b 


1.14 ± 0.08 


1.07 


B(L>* ^ L>°7r°) 


0.619 ± 0.029 


0.826 


Id*/ 




1.36 


B(L>* ^ D°7) 


0.381 ± 0.029 


0.174 


Ib* 1 f B 




1.10 


B(S -> i<r*7) 


(5.7 ± 3.3) 10 -5 


11.4 



Table 3: Calculated values of a range of observables not included in fitting the model's parameters. The 
"Obs." values are extracted from Refs. [15,17,18,24]. (Table adapted from Ref. [6].) 





Obs. 


Calc. 




Obs. 


Calc. 


Af^ D '(0) 




0.57 


AB^D* ( f B->D~\ 
■"■1 V-max ) 




0.88 






0.56 


AB^D' ( f B->D*\ 




1.16 


V B ^ D '(0) 




0.70 


v B ^ D '{t B z D ') 




1.47 




1.30 ± 0.39 


1.32 


R^ D '(i) 


0.64 ± 0.29 


1.04 






1.23 


R B ^ D '(w m!iX ) 




0.98 


a B ^P 




0.60 


fB^Dt f B^D\ 




1.21 


/f— (0) 




0.45 


fB^-K (j-B^1T\ 

J+ I, ''max I 




3.73 


A B ^ P (0) 




0.47 


\B^p,,B^p\ 
A \ I'-max ) 




0.45 


A B -"(0) 




0.50 


A B ^ p (t B ^) 




0.81 


V B ^»(0) 




0.68 


v B - p (C^) 




1.17 


i2f^(l) 




1.15 


R B b ^(1) 




0.80 


R B ^ p (w max ) 




1.44 


R B ^ p {w max ) 




1.06 



Table 4: Calculated values of some b — > c and b — > u transition form factor observables not included in fitting 
the model's parameters. The "Obs." values are extracted from Refs. [15,18]. (Table adapted from Ref. [6].) 



Takahashi-idcntity-preserving Bcthc-Salpeter ker- 
nel applicable to the study of heavy-meson masses. 
That would provide further insight into the struc- 
ture of heavy-meson bound state amplitudes, an 
integral part of these calculations for which only 
rudimentary models are currently available. It 
would also assist in constraining DSE phenome- 
nology via a comparison with calculations and 
models of the heavy-quark potential. 
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